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Abstract 

We derive a new, exact and transparent expansion for option smiles, which lends itself both 

to analytical approximation and to congenial numerical treatments. We show that the skew and 
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CN ■ the curvature of the smile can be computed as exotic options, for which the Hedged Monte Carlo 
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Q_i, method is particularly well suited. When applied to options on the S&P index, we find that the skew 
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and the curvature of the smile are very poorly reproduced by the standard Edgeworth (cumulant) 
expansion. Most notably, the relation between the skew and the skewness is inverted at small and 



q/ large vols, a feature that none of the models studied so far is able to reproduce. Furthermore, 

CL : 

the around-the-money curvature of the smile is found to be very small, in stark contrast with the 
t4-H ' highly kurtic nature of the returns. 
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Understanding the shape of volatility smiles in option markets is arguably one of the 



QQ. 



most active field of research in quantitative finance [1|, |2j. By definition, the existence 
of an option smile is the sign that the standard Black-Scholes model is not an adequate 
representation of the stochastic dynamics of financial assets. A huge variety of models 
have been proposed over the years in order to account for the non-Gaussian nature of 
price changes and the corresponding option smiles: jumps and Levy processes [3j, GARCH 
and stochastic volatility models 41 (including the popular Heston and SABR models |5j), 
multiscale (multifractal) models 6|, |7J, mixed jumps/stochastic vol. models, etc. Another 
model-free strand of research that sheds light on the origin and the general structure of 
option smiles, is to assume that the corrections to the Black-Scholes model can somehow 
be considered as "small" . A natural idea is to use a Edgeworth cumulant expansion of the 
distribution of the price change tt = T g ° (see Appendix). Working with the Bachelier 
model as benchmark (i.e. {S±)t>o follows a Brownian motion instead of a Geometric Brownian 
motion), the authors of 8|, |9| have obtained the following cumulant expansion for the smile: 

where ay/T denotes the standard deviation of rr, -M is the rescaled moneyness of the option, 
and St <C 1 and kt <C 1 are respectively the skewness and kurtosis of tt, which are assumed 



to be small for the expansion to make sense. 1 Note that the authors of 10J independently 
obtained a similar formula within a Black Scholes context, around the same time. The two 
formulas coincide for ayT <C 1, which is expected because in that limit the Black-Scholes 
model becomes identical to the Bachelier model. 

The interest of the above formula is that it allows one to understand why smiles have 
generically an asymmetric parabolic shape, with an asymmetry related to the skewness of the 
distribution of the underlying, while the curvature of the smile is proportional to its kurtosis. 
Since from general arguments the distribution of rv becomes Gaussian at large times (albeit 
perhaps very slowly due to long-memory effects [9j), the asymmetry and curvature of the 
smile are expected to go down with maturity T, as indeed observed in option markets. 

However, the expansion (TJQ) involves moments up to order 4 and as such is troublesome 
both theoretically and practically. Indeed, the cumulative distribution of returns -P(|rr| > x) 



1 Eq. ([I} in fact assumes that S|h <§; kt, which is often the case in practice. 



is found empirically to decrease as x~ At with // p» 3 for many assets (see e.g. Ill, |l2j) and 
thus the moment of order 3 is formally divergent. At any rate, it is in practice very difficult 
to estimate moments of order 3 and 4 because of the noise induced by large events. One 
possibility is to replace the skewness and kurtosis of tt by some lower moment approxima- 
tions (cf. |9(). However, this procedure is ambiguous as there is freedom in the choice of 
lower moments to be used. 

The purpose of this note is to establish a different smile expansion formula, which is 
general, rigorous and involves no moments of order greater than 2. This new smile formula 
lends itself to analytical treatments, which allow one to recover, for example, the recent 
results of Bergomi & Guyon [13| . But perha ps m ore importantly, our formula can be coupled 
to the "Hedged Monte-Carlo Method" of 14j, |l5[ to yield a powerful numerical method, 
which provides accurate estimates of the smile parameter for arbitrarily complex models of 
the underlying. One may even use historical data directly, short-circuiting any modelling 
assumption. Our formula, derived in the Appendix, reads: 



ass = cr(a T + PtM + ~f T M 2 + 0(M 3 )), 
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with the coefficients oct, Pt, It given by: 



(Xt 



\e[\u t \ 



-[1-2P(«t>0)], lT 



f^°)-2V ® 



where ut = rx/cry/T and pr(-) is the density of Ut- Using a Edgeworth expansion, one can 
show that both smile formulae ([1]) and ([2]) coincide in the limit where cumulants are small and 
cumulants higher than four can be neglected (see Appendix). One finds in particular flx ~ 
St/6. But what is remarkable is that while the "old" smile formula ([T]) is highly sensitive to 
extreme events, the new formula §Z§ only involves low moments of the distribution of ut- In 
particular, the skew of the smile, as measured by the coefficient /3t, is technically a moment 
of order 0, i.e. the large events do not play any role at all. 

This new smile formula can be used in different ways. One can for example obtain directly 
the coefficients a>T, @Ti1t i n a vol-of-vol expansion, recovering the Bergomi- Guyon results 
to lowest order 16j. One of the salient results of Bergomi & Guyon 13j is that to leading 
order in vol-of-vol and for a broad family of linear models, the Edgeworth expansion result 
(3t = St/Q holds, i.e. the skew of the smile j3, and the skewness St of the distribution of ut, 
are very simply related. A first exercise is to use our formula to investigate the case where 



the correlation between stock returns and volatility is non-linear. Assume for example a 
stochastic volatility model where Tt = f a t dW t) with dW t the standard Wiener process, 
and where the instantaneous volatility a\ is given by: 

a 2 = a 2 [1 + 2e(|£t + 0\ l €t+fl<0 )] , & = u [ e^^ dW u , e > 0. (4) 

Jo 

Here a 2 is a "baseline" volatility (which might itself evolve over long times scales, see below). 

When the threshold 9 is zero, the above equation means that when recent returns are 
ne gati ve, the future volatility will be higher than usual (this is the standard leverage effect 
[9|, |l7j), while when the recent returns are positive, there is no impact on the future vol. 
When 9 > 0, only relatively large negative returns will increase the volatility, while even 
small positive returns increase the volatility when 9 < 0. The coefficients (3t and St can be 
computed to first order in e in this case [la]. One finds that the equality (3t = St/6 only 
holds when 9 = 0. When 9 > 0, one the other hand, the skew amplitude \Pt\ is smaller than 
|iSt|/6, and vice-versa when 9 < 0. This shows that even for small vol-of-vols, non linear 
effects can significantly affect the relation between skew and skewness, and that there can 
be no general link between the two. 

Although analytical results are interesting, we believe that the true added-value of our 
new smile expansion comes from the following remark: the three coefficents «t, /3t an d 7r 
can be interpreted as the average payoff of some "exotic" options. Indeed, a^ is simply an 
at-the-money straddle, /3t is related to an at-the-money binary option, and the first term of 
-)t is a "no move" option which pays if the underlying ends very close to its initial price. 2 
Having interpreted these coefficients as option prices, one can use the "Hedged Monte-Carlo" 



(HMC) method proposed in 

nr 

at least two reasons (see [9|, [14 



14l 15j ] to price these options numerically. This is interesting for 
[l5|, llSI] for a more thorough discussion): first, the HMC has 



by construction a low variance, that allows one to price options accurately with a relatively 
small number of paths; second, the P&L of the hedge automatically transforms historical 
or empirical probabilities into risk-neutral ones. Technically, the pricing of the "exotic" 
options is done in a way very similar to what is described in [14j, ll5j, except that we do not 
try here to determine the optimal hedge for these options, but rather use a (sub-optimal) 



2 It is in fact common folklore, inherited from Eq. ([1}, that buying slightly out of the money and selling 
slightly in the money allows one to bet on the "skewness", while a butterfly trade allows one to bet on 
the "kurtosis". 



Black-Scholes hedge. This reduces HMC to a standard reduced-variance Monte-Carlo [Ji 
which is conceptually much simpler and faster numerically. Although only approximate, this 
pruned down version of HMC is accurate enough for the present purpose. The other subtle 
point concerns the pricing of the "no move" option, since one has to introduce a window 
with a small but finite width S, and carefully extrapolate the results to 5 = 0. We in fact use 



a Gaussian payoff function, exp(— u^/25 2 )/V2ii5 2 for that purpose. We have calibrated the 
method on known cases - for example on the above non-linear model, with good agreement 
between the analytical calculations and the HMC results. 

We now apply our method to real time series, using historical trajectories to price our 
exotic options and determine directly the skew 0? and curvature 77- of the smile, without any 
a priori model. The data we use is the time series of daily returns (open-high-low-close) of 
the S&P500 index in the period 1/1/1970 - 31/12/2011. We divide the sample into two bins: 
one corresponding to high volatilities (larger than the median), the other to low volatilities, 
where the volatility is a 20 day exponential moving average of a Rogers-Satchell estimator 
of the squared daily volatility before the day the smile coefficients are determined. We show 
in Fig. 1 the "fair" skew /3t as a function of T for options between 1 and 20 days, that we 
compare with the prediction St/Q based on a direct measure of the third moment skewness 
of the distribution of ut- We find, perhaps surprisingly, that while |<Sr|/6 is systematically 
larger than the skew |/3^| for large volatilities, the opposite is true for low volatilities. Let us 
emphasize here that we are not speaking of implied volatility smiles from option markets, 
but rather of theoretical predictions of what the fair parameters of the smile should be. 



Including some information from option markets could be done along the lines of 20]. 

In Fig. 2, we show the smile curvature 7^ as a function of T for the same two volatil- 
ity regimes above, and compare it to kt/24, the curvature obtained from the Edgeworth 
expansion ([!]). Here again the results are surprising: we find that the "fair" curvature 7^ 
of the smile around the money is close to zero, both for high and low volatilities. This is 
very different from what the empirical kurtosis of the distribution suggests, especially in the 
high vol regime where the kurtosis of returns is empirically very substantial, in agreement 
with the fact, noted above, that the distribution of returns exhibits power-law tails. 3 This 
is a striking illustration of the misleading character of the Edgeworth expansion ([T|): the 



The kurtosis might even be mathematically divergent. 
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FIG. 1: S&P500 data. Skew fix and skewness St/6 for different volatility regimes. Note the 
inversion of the order between the two quantities as the volatility increases. 
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FIG. 2: S&P500 data. Curvature 7t and kurtosis k^/24 for different volatility regimes. Note that, 
surprisingly, the curvature of the smile is found to be very small in both volatility regimes. 

around-the-money smile of index options should be close to a straight line with very little 
curvature, as indeed often seen on implied smiles, whereas the Edgeworth expansion would 
predict highly curved smiles. 

The results on the skew of the smile shown in Fig. 1 are clearly inconsistent with the 



relation (3t = St/6, and cannot be understood either within the above non-linear leverage 
effect with a fixed parameter 9. A way out would be to assume that 9 itself depends on 
the "baseline" volatility a in Eq. (TJJ above, with 9 > during high volatility periods and 
9 < during low volatility periods. This could be interpreted as meaning that at high vols, 
a significant down-trend is required to drive the volatility even higher, whereas at low vols, 
it takes a moderate down-trend to move the volatility up. It is also instructive to compare 
these results with the predictions of a popular asymmetric GARCH model used to model 



the dynamics of the S&P500 index [2l|. Writing the returs as a t r] t with rj t ~ A/"(0, 1), the 
dynamics of a t is postulated to be: 



' 2 1 



(5) 



a t = a(l + xt)] Xt+i = PXt + v{\ + Xt) 

We have run our HMC on synthetic time series generated using the asymmetric GARCH 
(GAARCH) model with parameters p = 0.9 and v = 0.1, corresponding to a "memory time" 
of the volatility equal to 10 days. One can of course now generate as many trajectories as 
needed to get good statistics. The results for /3t and 7t together with their corresponding 
cumulant expansion expressions are reported in Fig. 3 for the high volatility regime. We 
observe in this case a trend similar to the S&P500 data: the skewness |«St|/6 exceeds \/3t\ in 
absolute value and kt/24 is systematically larger than 77-, although the latter is significantly 
different from zero. However, for low volatilities (not shown) the behavior of the skew is 
qualitativaly different from the S&P500 data. In that case /3r and St/6 are very close 
(except for very small T), but with still |<St|/6 > \Pt\- In other words, the inversion 
observed on empirical data for low vols cannot be reproduced within the GAARCH model. 
That I St |/6 > \/3t\ can in fact be checked analytically in a small v expansion. Note this 



inequality implies that the "Skew Stickiness Ratio" introduced by Bergomi in [22j is expected 
to be smaller than 2 at small maturities. 

As a conclusion, we have set up a new, exact and transparent expansion for option 
smiles, which lends itself both to analytical approximation and, perhaps more importantly, 
to congenial numerical treatments. We have shown that the skew and the curvature of the 
smile can be computed as exotic options, for which the Hedged Monte Carlo method is 
particularly well suited. When applied to options on the S&P index, we have found that the 
skew and the curvature of the smile are very poorly reproduced by the standard Edgeworth 
(cumulant) expansion, or by the Bergomi-Guyon vol-of-vol expansion (at least to lowest 




FIG. 3: GAARCH synthetic data from equation fl5J), with p = 0.9, v = 0.1, and 10 6 Monte Carlo 
paths for high volatility regime. Upper panel: curvature 7t and kurtosis kt/24. Lower panel: 
skew /3t and skewness St/ 6- 

order). Most notably, the relation between the skew and the skewness is inverted at small 
and large vols, a feature that none of the models studied so far is able to reproduce. We 
have argued that some coupling between the leverage effect and the volatility level is needed 
to capture such non-trivial statistical features. Finally, the around-the-money curvature of 
the smile is found to be very small, in stark contrast with the highly kurtic nature of the 
returns. This would also require some more detailed understanding. We hope that the ideas 
and methods presented in this paper are of general interest, and provide a fruitful framework 
to interpret the information provided by option markets. 
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Technical Appendix 



Proof of the smile formula 



We denote the price St = S (l + r T ). We suppose that r T is centered and that a 2 T 
E[r T }. We set u T = r T /ay/T. We have: 

E[(S T - K)+] = S E[(r T -{%-- 1))+] = So f°° P(r T > x)dx 

b J K ~ S ° 

POO 

= SocrvT / P(u T > u)du 
Jm 



where we introduce the renormalized moneyness M = Jta ■ A direct expansion in M 
leads to: 

E[(S T - K) + ] « S aVT(^^ - MP(u T > 0) + ^pr(O)) 

In the Gaussian case, we get therefore: 

E[(S T - K)+] « S Q cr BS VT(4= " ^ + M } t- ) 
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where Mg — Tffcr is the Gaussian moneyness. 

We make now make the assumption: 



a BS = ff{p{ + (3'M G + i'M 2 G 



(6) 



so by definition M = Mq(oi' + 0' 'Mg + l'M G ). The smile formula corresponds to the 
following identity: 

, 1 



&BS\ 



Ma M% 1 



.^2 2 ^k' v 2 

Therefore, we get the following equation (to the order M G 

Mut\] 



. ( M. M P( l(r> 0)A(0)) 
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which finally lead to the following identification: 



a'M G +P'M G )P(u T > 0)4 
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Pr(0)). 
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/3 / = 7rE[|n T |](i-P(n r >0)) 
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<*' 2 pr(0) _ a' 
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+ nE[\u T \}a-P(u T >0)f 



Switching back to standard moneyness and using gbs = g{& + fiM. + 7.M 2 ) finally leads to: 



= % = V^(| - P(ut > 0)) 
7 = £-£ = V^(*# 
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The Edgeworth expansion 



Assuming St, kt <C 1, the Edgeworth expansion reads: 



-^ , 5 



P( Ji_ > x ) _ at(x) « -7V (3) (x) - — iV (4) (x). 



where iV(a?) = /_ 



-t^/2 



leads to the following approximations: 



uc ± 24 



dt, N = 1 — N and JV( n ) refers to the nth derivative of N. This 



\ - P(u T > 0) « -M TPr (0) « -^ < 1 

where M^ is the median of uy. With these assumptions, we get the following expressions 
for a, /?, 7: 

LI ~ ± 24 

/3« -Mr 

, ' 2 ^v^^ 1 + 8 I V2^ [ + 24^ ~ 24 

Within the aforementioned approximation, both smile formulas coincide since we know that 
the Edgeworth expansion also leads to: 



S T ~-6M T , «t«24(1- il-E[\uT\ 
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